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Abstract 

We develop a perturbation theory that describes bound states of solitons local- 
ized in a confined area. External forces and influence of inhomogeneities are taken 
into account as perturbations to exact solutions of the sine-Gordon equation. We 
have investigated two special cases of fluxon trapped by a microresistor and decay 
of a breather under dissipation. Also, we have carried out numerical simulations 
with dissipative sine-Gordon equation and made comparison with the McLaughlin- 
Scott theory. Significant distinction between the McLaughlin-Scott calculation for a 
breather decay and our numerical result indicates that the history dependence of the 
breather evolution can not be neglected even for small damping parameter. 







1 Introduction 



Solitons or solitary waves are ones of the most interesting objects in nature. Observation 
of a solitary wave on water was first documented more than one and a half century ago. 
Besides, solitons occur naturally in many other substances like optical fibre pQ, nonlinear 
lattices [2], hot and cold plasma [3] and are even mentioned responsible for Jupiter's red 
spots [3] and energy transfer in DNA [SJ. Most intensively solitons have been studied in 
long Josephson contacts. The matter is that switching from a superconducting to resistive 
state of the Josephson junction is related to appearance and motion of solitons in these 
contacts which are known also as Josephson vortices or fluxons. Such solitons or fluxons 
are well described by the sine-Gordon equation. 

In ideal case when the Josephson junction is infinitely long and narrow, Josephson soli- 
tons can be described analytically by well known exact solutions of the sine-Gordon equa- 
tion. However, there is always dissipation associated with quasiparticle current through the 
Josephson junction and inhomogeneities associated with its width and thickness. More- 
over, the real physical systems are always subjected to influence of external forces. All 
these factors may have significant impact on soliton behaviour. 

Although, the strictly one dimensional sine-Gordon equation is integrable [7J E], the 
perturbations to this equation associated with the external forces and inhomogeneities 
spoil its integrability and the equation can not be solved exactly. Nevertheless, if their 
influence is small, the solution can be found perturbatively. The perturbation theory for 
solitons was described in details by Keener and McLaughlin JU]- Later, in application 
to dynamics of vortices in Josephson contacts, the perturbation analysis of sine-Gordon 
equation has been developed by McLaughlin and Scott [TT] . 

In many applications there appears a need in localized oscillatory solutions of sine- 
Gordon equation. For instance, when Josephson vortex is pinned by an inhomogeneity 
or there is a bound state of vortex with an antivortex known as breather. Breather may 
appear as a result of collision of a fluxon with an antifluxon or even in the process of 
measurements of switching current characteristics |B] • The role of breathers is ambiguous. 
Depending on our expectations, they can be parasitic excitations or, vice versa, a good 
substance for generation of THz waves. Recently we have proposed a device that may 
deliberatly generate and trap breathers 0. 

There have been many theoretical and numerical studies dedicated to continuous sine- 
Gordon breathers |T21 EH El EE3 EH|- In particular, decay of a breather into fluxon and 
antifluxon induced by the external current was studied by many authors ^2J- Moreover, 
it was shown |T3] that a breather can be stabilized by ac drive even in presence of energy 
losses. Also, influence of the boundaries on a breather dynamics [14J has been investigated 
and quantization of its energy spectrum J3| has been predicted. 

Nevertheless, despite of numerous theoretical studies, the dynamics of a breather under 
dissipation has not been fully understood. McLaughlin-Scott theory gets overcomplicated 
when applied to non-trivial solutions such as breather, whereas its simplifications fail to pre- 
dict the correct dynamics. We have performed numerical simulations of breather dynamics 
and found that there is a significant discrepancy with the McLaughlin-Scott calculation. In 
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particular, it manifests itself in the dependence of the breather energy on time, FigO The 
thin line is the dependence following from McLaughlin-Scott calculation (formula (5.5) in 
Ref. JI]), the solid line represents our numerical simulations. This discrepancy stimulated 
us to look into this problem once again and develop perturbation theory that is designed 
especially for localized solutions of sine-Gordon equation. We have found that at the con- 
struction of such theory it is very important to take into account the history dependence 
of the breather evolution. Also, we have carried out direct numerical simulations with 
dissipative sine-Gordon equation. The numerical results appear in perfect agreement with 
our theory. 
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Figure 1: Dissipative dynamics of a sine-Gordon breather: Dependence of the energy H SG 
of a breather on time t calculated according to McLaughlin-Scott's formula (5.5) from the 
Ref. jTTj is presented by thin line. Dependence of the energy H SG of a breather on time t 
calculated by direct numerical simulations of sine-Gordon equation with damping is shown 
by the thick solid line. The damping constant is a = 0.01. 
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2 Perturbation theory for localized solutions 



The reason why McLaughlin-Scott formula for breather decay (the formula (5.5) taken 
from the Ref. fails to predict correctly the dissipative dynamics of a breather is the 
following. Breather is an oscillatory solution that is characterized by some "phase" that 
depends on history of the evolution. In their general formulation McLaughlin and Scott 
treat this difficulty introducing the history dependent term u{t') dt! and allowing ad- 
ditional time-dependent modulation of the free parameters (such as initial positions of 
fluxons or phases of breathers) in the non-perturbed solution fUJ 111) . The modulation 
of the free parameters is governed by additional differential equations. Obviously, this 
leads to additional complication because of the coupled differential equations for modu- 
lated parameters. Moreover, with such modulation the original solution no longer satisfies 
the non-perturbed sine-Gordon equation exactly and one needs also to modify the effective 
perturbation [H]. Here we describe a method that does not involve the modulation of free 
parameters, but correctly deals with time-dependent dynamics due to appropriately chosen 
anzatz of non-perturbed sine-Gordon solution. 

Consider solution of (l+l)-dimensional sine-Gordon equation 

<fxx + sin<f = (1) 

in the form ip(g(u)x, g(u)ut,u) with g{u) = 1/\A ± u 2 . Such parametrization is natu- 
ral for sine-Gordon solutions such as solitons and their bound states. The sine-Gordon 
Hamiltonian is a functional of the field variable ip, 

(2) 

Substitution of <p = tp(g{u) x, g{u) ut, u) gives the effective energy as a function of a single 
parameter u, 

H s e f f (u) = H SG [<p(g(u)x, g(u)ut,u)} 

The second argument of (p(g(u) x, g(u) ut, u) that we call here a phase T(t) = g(u) ut 
can be written in different ways, such as T{t) = g{u) udt! or T{t) = f* g{u)udt'. 
Obviously, in case of u independent on time these cases are equivalent and the choice does 
not make any difference. However, this definition of the phase is very important when 
taking into account the influence of perturbations, as will be shown below. 

In the presence of perturbations we assume that the dominant effect is to modulate the 
parameter u = u(t). In other words, with appropriate choice of u(t) we may satisfy the 
perturbed sine-Gordon equation 

<p' - <p xx + sin ip = ef 

by the function (p = ip(g(u(t)) x, T(t), u(t)). Here, we take the perturbation ef in a general 
form 

ef = - y^Hi S(x - Xi) siny? - 7 - a<p 



H 
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dx 
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In contrast to the case of constant u, the choice of the non-perturbative solution is not 
unique anymore. Indeed, depending on the choice of the phase T(t), we come out with 
different functions of t. We will show, that with the appropriate choice of the phase T(t) 
we may correctly describe the time evolution of localized sine-Gordon solutions in presence 
of perturbations. We describe the dynamics by a single modulated parameter u = u(t) 
without introducing additional modulation of the free parameters. This gives considerable 
simplification and improvement because the other free parameters such as initial location of 
solitons or initial phases of breathers remain fixed and do not result in auxiliary differential 
equations like those introduced in the Ref. [TT] . 
Consider the anzatz 

<p(g(u(t))x,T(t),u(t)) with T{t) = [ g(u(t')) u[t') dt' (3) 

Jo 

where function (p is an exact solution of non-perturbative sine-Gordon equation In 
further consideration we omit highlighting the explicit dependence of the functions u = 
u(t) and T = T(t) for typographical convenience. Obviously, the drawback of the time 
modulation of u affects the time derivative of (p, 

<P = ^j_(p(g(u) x,T,u) = (fig u ux + (f 2 gu + (p 3 u 

Where tpi, tp 2 and ip 3 are derivatives of if with respect to first, second and third argument 
correspondingly. As we consider localized solutions confined in some area \x\ < C, the 
term tpig'(u)ux is of the order 0(e). The third term also can be neglected as it does not 
contain explicit linear terms in x and t. Therefore, 

^ = V2 gu + 0(e) (4) 

that remains valid even in the limit of large times, t — > oo. Obviously, another choice 
of T(t) would spoil this equation with terms explicitly dependent on time t, e.g. for 
T(t) = g(u{t)) J Q u{t') dt' we would have 

(p(g(u)x,T,u) = (pig u iix + (p 2 gu + Lp 2 g u u I u(t') dt' + <p 3 u 

Jo 

that contain a non-zero term J Q * u(t') dt' proportional to t. Thus, in this case the dynamics 
would not be correctly described on large time scales, t — > oo. Mclaughlin and Scott 
overcome this problem introducing additional modulation of free parameters. 
Substituting (j3J) to (J2J) we obtain the effective energy as a function of u(t), 

^ G b(^(n(t))a:,T(t),n(t))] ~ H^(u(t)) (5) 

that is valid for any values of t. It is important to note that this expression coincides exactly 
with the effective energy of non-perturbed solution (JSJ) and depends on time inderectly only 
via u(t). 
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In the presence of external forces, we may write the full Hamiltonian, 

H = H SG + H p 

and take into account the dissipative perturbations affecting the energy dissipation rate JI] 

dH f°° 2 , 

cap dx 



dt 

The Hamiltonian H p serves to describe non-dissipative perturbations induced by external 
potential forces. This could be microshorts, microresistors or applied driving current. 



H 1 

Thus, 

dH SG 



- r (e 

J~oo \ i 



x — Xj) <f sin f + 70 + af 2 ] dx 



dt 

Substituting (JSJ), we obtain the equation for parameter u 

u ^ = — / ( y ] Hi 8(x — Xi) if sin ip + 70 + cap 2 ) dx (6) 



J-oo \ i 



where if = <f(g(u) x,T) and ip ~ ifT g{u) u should be substituted. This equation is coupled 
to the equation for T, 

f = g{u)u (7) 

In some cases it can be convenient to rewrite this system of differential equations for 
independent variable T, 



% = ~ I-oo (E< t* s (. x - x i)V sin <^ + 7 W + af 2 )dx 

§ = (g(u)u)- 1 



(8) 



The dynamics is described by u(T(t)), where T(t) is a inverse function of t(T). 

3 Pinning by a microresistor 

Let us consider the following ansatz for a single soliton solution 



f(g(u)) x, T) — 4 arctanexp (g(u)x — T) , g(u) = 1/y/l — u 2 
subjected to the attractive potential of a microresistor 

ef = —jji5{x) smf, ji < 



The energy of a soliton is equal to 



u 



From (jUJ) and (J7J) we obtain the next system of coupled differential equations, 

a = \n (1 - u 2 ) sech 2 (T(t)) tanh(T(t)) 
T 



(9) 



We have found that after some simplifications, the McLaughlin-Scott's formula (4.3) from 
Ref. [TTj can be reduced to the exactly the same system of differential equations. Although, 
both approaches lead to exactly the same result, McLaughlin-Scott's formulation is, obvi- 
ously, more cumbersome. 
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Figure 2: The dependence of the speed u(T) on the phase T calculated using perturbation 
theory for localized sine-Gordon solutions at different damping rates. The top curve cor- 
responds to the damping constant is a = 0.05, the middle curve - a = 0.01 and the lower 
curve a = 0.005 . 
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Figure 3: The dependence of the time t on the phase T calculated using perturbation theory 
for localized sine-Gordon solutions at different damping rates. The top curve corresponds 
to the damping constant is a = 0.005, the middle curve - a — 0.01 and the lower curve 
corresponds to the highest damping a = 0.05 . 



4 Decay of a breather 

Consider a breather solution 



<f(g(u) x, T,u) — 4 arctan 



sinT 



u cosh (g(u)x) 

with g{u) = l/y/i + u 2 . As a perturbation we consider the dissipative term 

ef = -a(p 

The effective energy is 
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Figure 4: Dissipative dynamics of a sine-Gordon breather: Dependence of the energy 
H SG of a breather on time t calculated using perturbation theory for localized sine-Gordon 
solutions (thin line). Dependence of the energy H of a breather on time t calculated by 
direct numerical simulations of sine-Gordon equation with damping is shown by the thick 
solid line. The damping constant is a = 0.01. 



From (jHJ) we obtain the next system of coupled differential equations 



da _ (1+u 2 ) 3 / 2 eos 2 T 
dT a sin 2 T+u 2 



1 + 



u 2 arctanh( . si " T = 

\ Vsin 2 T+u 2 



sinT ysir 



(10) 



dt 
dT 



where u = u(T(t)). This is a new result that may not be obtained from the McLaughlin- 
Scott theory by straightforward manipulation. The system can be solved numerically. For 
an illustration we present the solution of these equations on Figs EJ and H3 There the 
initial conditions are taken as w(0.1) = 0.1 and £(0.1) = 0.1. One may notice the steplike 
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character of the dependence u(T), see Fig|21 The size of the steps increases with dampling 
indicating importance of the introduction of the phase T(t). This phase has also nontrivial 
dependence on time t. Its inverse function t(T) is presented in FigCJl One may notice that 
at some values of T there is a fast change of the slope. Obviously this is related to the 
steplike character of the dependence u(T). The dissipative dynamics of a breather is also 
well reflected by time dependence of its energy, Fig0] The results are in perfect agreement 
with our numerical simulations using the complete sine-Gordon equation with dissipative 
term. 

5 Conclusion 

In summary, we have found that our perturbation theory describes well the dynamics 
of localized excitations subjected to influence of external forces such as various inhomo- 
geneities and damping associated with quasiparticle current. In particular, we have de- 
scribed a fluxon trapped in a potential well which could be related to a microresistor in 
the Josephson junction. Here the equations derived with the use of our method coincides 
identically with equations derived by McLaughlin-Scott However the derivation of 

these equations obtained by our method is significantly simpler. Second, we have described 
the decay of the breather under dissipation. In this case, the equations are different from 
the McLaughlin-Scott's [O]. According to our calculation, the breather is decaying signif- 
icantly faster. In order to resolve this difference we have performed numerical simulations 
with dissipative sine-Gordon equation. The results of these numerical simulations are in 
perfect agreement with our theoretical results. The comparison of our perturbation theory 
with the McLaughlin-Scott's calculation [llj indicates that the history dependence of the 
breather evolution has a strong influence on its dynamics even at low damping. 

To conclude, we have developed the perturbation theory which perfectly describes the 
localized in space solutions of sine-Gordon equation. Such study can be important for new 
devices such as fluxon collider 5| or other devices based on newly discovered flux cloning 
effect ^7] where the dissipative dynamics and the breather excitations may play a key role. 
The theory may allow generalization to higher dimensions. This can be of use to study 
localized pulsating solutions of sine-Gordon in two spacial dimensions [TH] . 
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